Buckling induced Zener polaron instability in half-doped manganites 
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By calculating the electronic, orbital, magnetic and lattice structure of half-doped manganites we 
establish the central role of oxygen buckling caused by the tilting of the MnOe octahedra, which 
is, in essence, a steric effect. The buckling itself does not change the Mn-Mn bonds. Instead, 
it drives the system toward an instability where these bonds can dimerize. In the presence of 
electron-electron interactions, this instability can fully develop and beyond a critical buckling a 
Zener polaron groundstate with dimerized spins, lattice, and orbitals forms spontaneously, resulting 
in a multiferroic state. 

PACS numbers: 75.47.Lx,71.10.Fd, 75.85.-|-t, 75.25.Dk 



Materials where several degrees of freedom, spin, 
charge, orbital and lattice, are simultaneously active dis- 
play very rich phenomenologies and exotic ground states, 
including ferromagnetic metals or antiferromagnetic in- 
sulators with complex charge and orbital ordered states^. 
Manganites clearly belong to this class, as is evident from 
the observed colossal response to external fields of mag- 
netic, electric or strain originSri,^ which indeed reveals a 
subtle balance between different interactions. Despite the 
enormous experimental activity and the considerable the- 
oretical progress in understanding their rich phenomenol- 
ogy, manganites continue to unveil unexpected and fas- 
cinating properties like the coexistence of spontaneous 
ferroelectric polarization and magnetic order'^^^, which 
defines them as potential multiferroics, offering new op- 
portunities for, e.g., multifunctional magnetoelectric de- 
vicesSii^. 

In particular doped rare-earth manganites, 
i?i_a;Aa;Mn03 (i? and A are rare-earth and alkaline- 
earth cations), were recently suggested to have a 
potentially large multiferroic coupling—, particularly 
the half-doped Pri/zCai/zMnOa (PCMO)iiii^ and 
Lai/2Cai/2Mn03 (LCMO)i^. Apart from the possible 
appearance of a multiferroic state, the physical proper- 
ties in the low-temperature regime of these compounds 
are still debated. It is widely accepted that the local Mn 
magnetic moments can order in the so-called CE-type 
antiferromagnetic (AFM-CE) configuration character- 
ized by double-zigzag spin chainsi^ (Tc = 175/155 K 
for PCMO/LCMO), but the picture is less clear for the 
charge (CO) and orbital ordering (00) that are observed 
and two possible scenarios have been envisaged. On the 
one hand a checkerboard pattern (CB) of alternating 
Mn'*+ and orbital-ordered Mn'^"'" has been proposed in a 
low-symmetry P2i/m space group^^. On the other hand, 
a Zener polaron picture (ZP) has been suggested in the 
low-symmetry P2inm space group, where equivalent 
Mn^"*" ions couple into ferromagnetic dimers, displaying 
an 00 pattern with orbital states aligned along the 
same direction within each dimer. In this picture the 
CO originates from holes on oxygen atoms that bridge 



the Mn dimersiiii^. The symmetry of the ZP state 
with P2inm structure allows for a spontaneous electric 
polarizatiortii, which is related microscopically to the 
presence of inequivalent oxygen sites^ii^. 

From the theoretical point of view, the CB picture has 
received considerable attention in recent years and has 
been investigated by means of computational and mean- 
field analysis of realistic models^. The AFM-CE con- 
figuration was shown to give rise to a band-insulating 
behavior that emerges from the topology of the zigzag 
chains in the presence of large Hund coupling between 
t2g Mn core spins and Cg electron o^'^^i^^ , which is also re- 
sponsible for the emerging orbital-ordering pattern. The 
checkerboard CO pattern can then come about via ei- 
ther electron-electron'^'' or electron-lattice Jahn- Teller- 
interactions. On the other hand, the ZP state has been 
studied mainly in the framework of first-principles calcu- 
lationsi2iHiiSi2£ due to its complex nature and because of 
the expected importance of the oxygens ions, which are 
usually not included in model Hamiltonian approaches to 
manganites. 

In this Brief Report we consider a straightforward 
extension of the degenerate double-exchange (DDEX) 
model that includes the oxygen-buckling mode modulat- 
ing the Mn-Mn hopping processes by changing the Mn-0- 
Mn angle (p. Such a buckling is ubiquitous in perovskites 
and is induced by a tilting and rotation of the MnOg octa- 
hedra. The amount of tilting depends on the ionic radius 
of the rare-earth cation that is used and is, in essence, a 
steric, and thus very robust, effect. We find that the in- 
terplay of this buckling mode with the electron-electron 
(e-e) interaction generates an instability toward Mn-Mn 
dimerization. The emerging dimerized state has the qual- 
itative features predicted for the ZP state, in particu- 
lar, the same 00 pattern and a strongly reduced charge 
disproportion between Mn ions. The inclusion of fur- 
ther electron-lattice (e-l) interactions, in particular, the 
coupling with the Jahn- Teller (JT) mode, further stabi- 
lizes this phase. We analyze the different contributions 
to the ferroelectric polarization of the multiferroic state 
that emerges. 
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Model. In the infinite Hund coupling approximation, 
Cg electrons effectively only hop between sites with fer- 
romagnetically aligned core spins. As a consequence in 
the AFM-CE phase hopping is allowed merely within one- 
dimensional zigzag chains of parallel core spins, changing 
direction as {..x, x, y, y..}. In the model Hamiltonian we 
also include the effect of the GdFeOs-type distortions, a 
change of the Mn-O-Mn angle if, on the hopping param- 
eters as well as the e-e and e-l interactions: 

^ = - E ^% ^-'^w + f/ E + (1) 

i ini 

The first term corresponds to the oxygen-mediated elec- 
tron transfer between nearest-neighbor manganese sites 
along the {...x,x,y,y...} chain, where dj ^ creates a 
particle at site i in electronic states stemming from 
eg orbitals d,j.-2_y2 [a ~ a) and d^z^^r^ (a ~ b) and 
taa = Stocos^ 'fitj/4:,tll = toCOSLpij/A and f^^^ = 
±V3to cos^ (pij / 4^, where to = {pda^ is taken as the 
bare hopping parameter. The Mn-O-Mn angle depends 
on the oxygen positions as coi{ipi±i/2) = luij^^lii where 
^^i±i/2 is the displacement of the oxygen atoms around 
the ith Mn site perpendicular to the bond, see Fig. HJa). 
The sign in front of the interorbital amplitude is 
determined by the direction of the vector connecting 
two neighboring sites along the zigzag chain; we can 
distinguish then between bridge sites, where hopping 
does not change direction, and corner sites, where it 
does. The second term describes the Coulomb inter- 
action, and the last terms describe the e-l interaction, 
with coupling A and operators ni — dl^d^^ + dl^df^^, 

Ti = dlidai - 4>Ai and rf = d^d^^ -h 4^c^a^■ <9™ 
represents the different phonon modes, i.e. breathing, 
Jahn- Teller and buckling modes, with stiffness constants 
Km- The buckling mode on the chain corresponds to 
Qi.buckle = {ut+i/2 +'"j-i/2 - (-l)*2uo), where uq is 
the bare value of the buckling distortion. We adopt 
standard values for the stiffness constant of the breath- 
ing and JT modes, ki = 2, K2 = K3 = kjt = 1 and 
use Kbuckic = 10-- In the following we express the 
strength of the e-l coupling in units of the static JT en- 
ergy EjT = AV(2/«jt<o)- 

Calculations. We determine the ground state wave- 
function variationally with respect to the phonon modes, 
treating the on-site interaction at the mean-field level. 
Since the Coulomb term appearing in Eq. ([1]) is ro- 
tationally invariant in the orbital space, the decou- 
pling Hartree-Fock procedure is implemented in the lo- 
cal optimum Cg-electron basis, which is determined self- 
consistently. This allows us to retain the orbital de- 
grees of freedom in the mean-field framework; should 
we linearize the interaction term in the original ba- 
sis, we would fix the orbitals to be either x^ — y^ or 
3_j,2 _ j,2 ^ -jij^g linearized interaction term is 4 riia riib ~ 
2 in,{n,) - Tzi{Tzi) - T^i{Tzi) ) - {n^f -\- {Tzi)"^ + {T^i^ . We 
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FIG. 1: (a) Schematic representation of the buckling distor- 
tion in the one-dimensional zigzag chain and of the dimeriza- 
tion eflFect induced by the local interaction (manganeses with 
opposite spins belong to different chains), (b) Energy gain of 
the dimerized (ZP) state with respect to the non-dimerized 
(CE) state as a function of U (at Ejt = 0) for an initial 
distortion uq — 0.3. Inset: equilibrium distortion u^q for the 
non-dimerized state compared to the staggered one found in 
the ZP state, (c) Difference of the on-site orbital occupancy 
between the most and least occupied orbitals as a function of 
U in the dimerized state, on corner (circle) and bridge (tri- 
angle) sites. Inset: charge disproportion between corner and 
bridge Mn sites. Dotted lines represent the same quantities 
in the non-dimerized state, (d) Schematic representation of 
the orbital-order pattern in the dimerized state. 



first consider the A = {Ejt = 0) limit. When the ini- 
tial buckling is small {uq < 0.2), we recover the standard 
resuh for the CE-type DDEX modelAiZ. At half-fiUing 
the ground state of the non-interacting system is a band 
insulator with bridge Mn ions displaying (partially filled) 
orbitally ordered states 3a;^— r^/3y^— r^, while corner Mn 
ions do not display any 00. As mentioned above, such 
a band insulating state has a topological origin which is 
ultimately related to the constraint on eg electron mobil- 
ity introduced by the underlying magnetic configuration 
of t2g spin o^^d^ . Since both orbitals of the local electron- 
state set are partially occupied at corner sites but only 
one is occupied at bridge sites, the inter-orbital Coulomb 
interaction U acts differently on them, inducing a charge 
transfer from more correlated (corner) to less correlated 
(bridge) sites, i.e. U induces a CB charge-order pattern 
with Mn+'^-^~''/Mn+^-^"'"* on bridge/corner sites. As a 
consequence of the competition of hopping and elastic 
energies - the kinetic energy is of course minimal for the 
undistorted structure {(p = 180°) - the optimal oxygen 
displacements along the chains are Ui±i/2 — u^q < uq. 

For larger initial distortions, however, the system gains 
energy for physical values of U by dimerizing, i.e. the 
optimal values of Uj^x/2 become different at even and 
odd bonds causing an effective hopping dimerization, as 
shown in FigsHI^a) and[ljb). We refer to this dimerized 
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FIG. 2: (left) Evolution of total energy as a function of ii_ — 
{ui — ii2)/2 at fixed optimum u+ — (ui + U2)/2 at [/ = 
and U — 1 (EjT = 0); the inset shows the onset of the two 
contributions to polarization at U = 1. (right) Polarization 
evaluated at uq — 0.3 and Ejt = as a function of U (the 
electron charge e is set to unity). 
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FIG. 3: (left) Difference in on-site orbital occupancy as a 
function of Ejt in the dimerized state, on corner (circle) 
and bridge (triangle) sites. Inset: charge disproportion be- 
tween corner and bridge sites. Dotted lines label the non- 
dimerized results. Open symbols show the same quantities 
for the breathing mode switched off. (right) Calculated po- 
larization as a function of Ejt at uo = 0.3, U — (solid and 
open symbols: with and without the breathing mode). In- 
set: energy gain of the dimerized phase with (solid line) and 
without (dashed line) breathing mode. 



phase as a Zener polaron state, even if we are at this point 
not considering the magnetic moments of such dimers. 
Correspondingly, also an orbital ordering establishes at 
corner sites, with the most occupied orbital aligned along 
the direction in which the hopping amplitude is larger. 
A good measure of the orbital ordering is Ago = Pa — Pi, 
where pa — (il't'^a) the local average occupancy of the 
orbital state a and 4'a — -~ sm{9a/2) da + cos{9a/2) di,. 
Here a, b label the local most and least occupied or- 
bital states, which we find to be — (3y^ — r^) 
and 2/^ — z"^ {x^ — z^) respectively. As a secondary ef- 
fect, the correlation- induced charge transfer from corner 
to bridge sites is considerably reduced with respect to 
the non-dimerized state (see inset of FigHJc)), a direct 
consequence of the small energy cost for an electron to 
partially fill a single orbital state on corner sites. This 
suggests that the gain in electronic energy comes from 
the minimization of the on-site Coulomb energy at corner 
sites, directly related to the onset of the orbital ordering. 
The 00/CO patterns which characterize the dimerized 
state in the proposed model agree perfectly well with that 
expected in the proposed ZP picture for the half-doped 
manganites. We note also that the Mn-O-Mn angle de- 
pends on the nature of the bond (parallel vs. perpendicu- 
lar orbital states along each chain, ferromagnetic vs. anti- 
ferromagnetic between different chains); in the dimerized 
phase we end up with three different angles, correspond- 
ing to the involved bonds, in qualitative agreement with 
the PCMO measured anglesii. 

Multiferroicity. Having established that buckling of 
the MnOg octahedra drives the system toward an insta- 
bility where the dimerized Zener polaron state can form 
in the presence of electron-electron in interactions, we 
now focus on the possible ferroelectricity of the resulting 
ground state. The bridging oxygens located in the middle 
of the bonds are no longer inversion-symmetry centers in 
the dimerized state, allowing for an in-plane ferroelectric 
polarization Pion perpendicular to the chains direction 
(see FiglTJa)). In analogy with the displacement mecha- 



nism induced by the exchange coupling in undoped rare- 
earth manganites, we refer to the mechanism that dif- 
ferentiates the bridging oxygen positions in the proposed 
model as a correlation-induced orbital-striction. One can 
estimate Pion as the oxygen displacement per "unit cell" 
containing one Mn and two O, averaged over the true 
unit cell along the chain, consisting of four Mn sites. But 
in addition a purely electronic contribution to the polar- 
ization is generated by the interplay of the orbital order 
and the specific topology of the chains. In fact, eg Bloch 
electrons acquire a Berry phase in the 00 state, and a 
different geometric phase coming from the change of the 
hopping directio n^^i^^ , whose interplay gives rise to a fi- 
nite polarization^^. In the non-dimerized AFM-CE state, 
the phase associated to 00 appears only on bridge Mn, 
whereas the hopping phase due to the zigzag geometry of 
the chains originates at corner Mn, thus giving no polar- 
ization—. On the other hand the onset of 00 on corner 
Mn sites allows for a finite P, whose magnitude is directly 
related to the strength of the orbital order. We evaluated 
this Berry phase contribution to the polarizatio n^'^'^"^ , 
finding a sizable polarization parallel to the chains di- 
rection, in excellent qualitative agreement with ab-initio 
calculations for both LCMO and FCMOi^^H. Both the 
electronic and ionic contributions for the one-dimensional 
chain are shown in Fig 121 

The same analysis can be repeated in the presence of 
the e-l interaction. On the basis of our previous find- 
ings, we expect that the JT coupling, that may induce 
the occupation of a preferential local orbital state, would 
favour the onset of the ZP dimerized phase, whereas the 
breathing mode would probably suppress the tendency 
to dimcrizc, favouring instead a CB charge-ordering pat- 
tern. Indeed, at J7 = we find a very narrow win- 
dow of values for Ejt in which the dimerized phase sets 
in, with a very small energy gain with respect to the 
non-dimerized state. This energy gain comes from the 
JT coupling, which is thus less effective than correlation 
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FIG. 4: Phase diagram in the Ejt-U space at (left) uo = 0.3 
and (right) uo = 0.25; ZP is the dimerized state and CB the 
non-dimerized one. 



in inducing the differentiation of oxygens displacement, 
even if the resulting picture is qualitatively the same (see 
Figs 121 left, and 121 right). In fact, when the breathing 
mode coupling is switched off the ZP stability region is 
significantly enlarged, a further indirect proof that the 
orbital ordering on corner sites plays a relevant role in 
the onset of the dimerized state. On the other hand, 
the e-e and the JT e-l coupling cooperate to enhance the 
orbital-striction mechanism, even in the presence of the 
breathing mode, as clearly shown in the phase diagram 
in the Ejt-U space shown in Fig HI 

Conclusions. The instability of the system toward 
the formation of a ZP dimerized phase depends crucially 
on the competition and subtle balance between double- 
exchange, electron correlation and the lattice elastic en- 
ergy. Repeating the analysis above for different values 
for fcbuckic, we find that the dimerization becomes weaker 
for smaller fcbuckic, and is thus more stable for a stifFer 



lattice. Indeed, in the first case the hopping energy gain 
due to the increase of (p at both even and odd bonds 
overwhelms the energy gain that would result from these 
bonds becoming inequivalent; on the other hand, the 
orbital-striction mechanism comes into play when the 
Mn-O-Mn angle cannot be increased very much at any 
bond. The initial degree of buckling is essential: the more 
distorted the system is, the more effective the Coulomb 
and JT interactions are in stabilizing ZP, see Figl4l In 
our model the magnetic degrees of freedom of the dimers 
and the charge degrees of freedom of the oxygens are 
not included, which is justified by the observation that 
the most stable magnetic configuration for the half-doped 
case is very close to the AFM-CE state^^, even if the ge- 
ometry of the exchange couplings between dimers, that 
is of triangular rather than square type, can stabilize 
magnetic configurations different form AFM-CE thus 
lessening the instability^^. The oxygen charge degrees 
of freedom can become relevant because the enhanced 
transfer energy inside each dimer can activate a charge 
transfer from bridging oxygens to manganeses, which is 
expected to compensate for the charge disproportion be- 
tween the corner and bridge Mn sites in our model. This 
can change the value of the ferroelectric polarization but 
will not affect our conclusion on the presence of a Zener 
polaron instability that is induced by the buckling and 
rotation of MnOg octahedra. The qualitative features of 
the emerging ZP state are reproduced very well in our 
model Hamiltonian for manganeses only; in the present 
context a charge transfer from oxygen to manganese is 
rather the consequence of the instability. 
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